BNL-HET-98/12, TTP98-14, |hep-ph/ 9804215 

Two-loop QCD corrections to semileptonic b decays at an 

intermediate recoil 



Andrzej Czarnecki 
Physics Department, Brookhaven National Laboratory, 
Upton, NY 11973 



Kirill Melnikov 

Institut fur Theoretische Teilchenphysik, Universitat Karlsruhe, 
D-76128 Karlsruhe, Germany 
(March 1998) 



Abstract 



We present complete 0{a 2 s ) corrections to the decay b — ► clvi at the point 
where the invariant mass of the leptons y/q 2 equals the c quark mass. We 
use this result, together with previously obtained corrections at the ends of 
the q 2 spectrum, to estimate the total width of the semileptonic b — > c decay 
with 0(a 2 ) accuracy, essential for the \V c b\ determination. 

Precise determination of the Cabibbo-Kobayashi-Maskawa (CKM) matrix parameters 
is one of the main goals of the present and upcoming high energy experiments. One of 
the three parameters which can be measured with the highest accuracy is \Vcb\, for which 
there are two complementary methods (for a recent review, see e.g. fHH). The first one is 
based on studying the exclusive decay B — > D* + leptons near the zero recoil end of the 
phase space. The other, so-called inclusive method examines the total semileptonic width 
of the B meson, T s i(i? — > XJui). For both determinations the experimental precision has 
reached a few per cent level and will be further improved in the future I?-factories. To fully 
exploit those measurements, the theoretical prediction for the b quark semileptonic width 
must be known with comparable precision. This requires a study of both short- distance 
(perturbative) and long-distance (non-perturbative) QCD effects. For the exclusive method 
the 0(a 2 ) corrections were calculated in ||. Technically, the most challenging part of that 
calculation was the QCD correction to the axial b — > c current normalization at zero recoil, 
obtained in |4|U and confirmed in a recent study || . Theoretical accuracy of the exclusive 
method is limited by the errors in the non-perturbative matrix elements which are enhanced 
due to the not so large mass of the c quark. 

In the inclusive method the non-perturbative corrections are somewhat smaller, with 
suppression by 1/m 2 . They are estimated |7[ to decrease the meson semileptonic width 
by approximately 5% compared to the free quark decay rate. With the non-perturbative 
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corrections under control, the perturbative corrections should be carefully analyzed. Ideally 
one would like to know the two loop QCD correction for an arbitrary invariant mass squared 
of the leptons (q 2 ) in the decay b — > clvi. This has been achieved for only a subset of correc- 
tions, enhanced by a large factor describing the running of the strong coupling constant |§ ; 
these so-called Brodsky-Lepage-Mackenzie (BLM) corrections || have been even resummed 
to all orders [|IJ. However, in absence of a complete calculation of the two-loop QCD effects, 
the related uncertainty in \ V c b\ determined from the inclusive method had to be guessed and 
various estimates were given in the literature. It is clear that the problem can be only solved 
by an explicit calculation of the complete 0(a 2 ) effects and therefore the need for such a 
calculation was often emphasized. In the present paper we demonstrate that, with the new 
result presented here, perturbative effects can be estimated with sufficient accuracy. 

For an arbitrary value of q 2 (and, by the same token, for the total width), the calculation 
of the second order corrections is a very difficult technical task. In earlier works we presented 
such results for the two extreme points, q 2 = JTTI] and q 2 = {m b — m c ) 2 Hfl. While at 



both points the non-BLM corrections were below 1% of the total semileptonic width, it was 
not clear what happens at intermediate values of q 2 , where a complete calculation appeared 
infeasible. Nevertheless, from the results in zero and maximal recoil points it was conjectured 
in [IT|, that the second order non-BLM corrections to the semileptonic decay width b — > clvi 
do not exceed 1%. 

Since then, we have found another kinematical point, q 2 = m 2 , where the complete 
calculation of the second order QCD corrections is possible. It is fortunate that for the 
physical value of m c this roughly corresponds to the middle of the q range, so that together 
with the endpoints we have nicely distributed information (see fig. [I]). 

Thus, the purpose of this paper is twofold. We first describe a calculation of the 0(a 2 ) 
correction to the differential semileptonic width of the b quark for q 2 = m 2 . We then use 
the results obtained in this paper, as well as in to estimate the 0(a 2 ) correction to 

the semileptonic decay width of the 6-quark. 

We start with a short description of the new calculation of the 0(a 2 ) correction for 
q 2 = m 2 . A complete description of the technical details of the present calculation will be 
presented elsewhere. The basic idea is the expansion of the decay amplitude in the velocity 
of the massive quark in the final state. 

The Feynman diagrams which describe this process can be divided up into three classes 
according to the number of real gluons emitted. In the first group there are the purely 
virtual two-loop corrections. Their value is known at the point m c = q = mt ) /2 from our 
previous study on the zero recoil line ||. Since the actual value of m c is not equal m b /2, we 
expand the diagrams in the parameter /3, related to the difference of these two quantities 
(see below) . The coefficients of this expansion are given by diagrams with higher powers of 
propagators. By solving a system of recurrence relations, obtained by integration by parts, 
we found an algorithm with which all these Feynman integrals can be expressed in terms of 
four non-trivial functions, given below in eq. (|J), and some simpler (single-scale) Feynman 
integrals. This is the most challenging part of this calculation and it limits the number of 
terms we have been able to obtain using the present computing facilities. 

In the other two groups there are diagrams with two real gluons or one real gluon and 
a virtual loop. It is relatively easy to perform calculations in these cases using (slightly 
modified) computational techniques described in detail in [13]. This possibility is related to 
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the fact that for q 2 = m 2 the three- and four-particle phase space integrals can be expressed 
in terms of the Euler T function. 

We write the differential semileptonic decay width of the decay b — > civ at q 2 = m 2 as 
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where ABorn,i,2 describe the m c /mb dependence in various orders in the strong coupling 
constant. Both Aeorn = (- 



1 — m 2 /m 2 ) J 1 — 4m 2 /m 2 and Ai are known in a closed analytical 
form [13,14 . A2 is calculated in the present paper. For the purpose of presentation we divide 
it up into four contributions according to the color factors: 
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The last term, Ah, describes the contributions of the massive b and c quark loops. Top 
quark contribution is suppressed by a factor ~ m\jm\ 
SU(3) group the color factors are C A = 3, C F = 4/3, Tr = 1/2. N F 
the quark flavors whose masses have been neglected (u, d, and s). 



m\jm 2 t and has been neglected. For the 

3 is the number of 



We have computed 5 terms of the expansion in (3 = (1 — 4m 2 /m 2 ) of A 
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In order 

to save space, in the formulas presented here we give only numerical values of the coefficients 
of the second and higher powers of (3 (these coefficients contain ln(/3); we use (3 = 0.64 to 
evaluate them). We use the MS scheme for the strong coupling constant and normalize 
it at the scale ^mbm c , which has been proposed as optimal for the semileptonic b — ► c 
transitions, in view of the proximity of the small velocity limit |[5 . Using the pole mass of 
the b and c quarks and expressing the one- loop corrections in terms of a s ( v /m{,m c ) we find 
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In the above formulas we use the notations /13 4 and R 2 for the values at uj = 1/2 of so 
denoted functions, for which analytical formulas are given in |J. Numerically they give 

fx ss 3.24460, / 3 « 12.3201, 

/ 4 « 7.83195, i? 2 « -0.72946. (4) 

We now would like to estimate the uncalculated remainder of the series given in @ 
and the error in the final result. In a series dnft n the remainder is less than the last 
known term multiplied by ft/(l — ft), provided that a n is a decreasing sequence, which we 
assume here. Therefore we estimate the final result by adding to the known terms half of 
ft I (1 — ft) times the last term. This additional term also gives a conservative estimate of the 
error. This procedure does not apply directly to Ap, where the coefficients appear to grow. 
From previous experience with similar calculations we think that this is because the series 
describing A^ has not yet achieved its asymptotic behavior ~ 1/n. For example, in the 
maximal recoil case [T2|] the first five terms of the expansion for behave rather wildly; 



nevertheless, they approximate the final result with an accuracy of about 25% (fortunately, 
the final result for A2 is rather insensitive to this error). We assign a similar error bar to 
A. jr. For ft = 0.64 we find: 

A A = -2.37(5), Ajr = 1.78(50), 

A L = 1.09(4), A H = -0.29(1). (5) 

Finally, we get for the correction defined in eq. @ (we add the errors in quadrature) 

A 2 = -4.72(14). (6) 

We now summarize the information about perturbative corrections to b —>■ c transitions and 
estimate the 0(a 2 s ) correction to the total semileptonic decay width of the 6-quark. We use 
the results presented in this paper, as well as in our previous papers P],[5|, HJ] . For arbitrary 
q 2 , we define: 
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where my, refers to the pole mass of the b quark. 

The BLM part of A2 was obtained in ||. The difference between the complete A2 and 
^blm - g ca n ec j a non~BLM correction. We introduce a quantity £, 

A 2 (q 2 ) - A^(q 2 ) 

aq ] ~ ' (8) 
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which describes the size of the non-BLM correction relative to the Born term as a function 
of q 2 . The values of £(q 2 ), for q 2 = fTlH , q 2 = m 2 (this paper) and q 2 = g^ ax = {m b — m c ) 2 
HH are, respectively, 0.65, 1.0, 0.06 (for m c /m b = 0.3). 

These three numbers permit us to estimate the 0(a 2 ) correction to the semileptonic 
decay width. For this purpose, we interpolate the non-BLM corrections with the function 
£(<? 2 ) — a 2<? 4 + &\q 2 + 0,0, where the coefficients a« are determined from the above values of 
£(0), £(m 2 ), £(<?max)- The function ^(g 2 ) is then integrated over q 2 , using the Born differential 
rate as the weight. As the result one gets an estimate of the non-BLM correction for the 
total semileptonic decay rate: 
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The validity of this estimate can be checked by using the same procedure to obtain 
an estimate of the BLM correction for the total decay rate using the known results at 
q 2 = 0, m 2 , (rrib — m c ) 2 as an input. 

Such a fit results in the value of the BLM correction — 8.6(a s /7r) 2 , to be compared with 
the exact result || — 9.8(a s /7r) 2 .[] We conclude that the accuracy of our simple fit is not 
worse than 30%, which also includes the errors in the three input data points. We therefore 
obtain the following formula for the total semileptonic decay rate r s i(6 — > cli>{): 
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where F(x) = 1 — 8x — 12x 2 ln(x) + 8x 3 — x 4 and we used m c /mb = 0.3. For the sake 
of clarity we separated the BLM and the non-BLM parts of the second order corrections. 
We also explicitly indicated the uncertainty in our estimate of the second order non-BLM 
correction. 

In principle, Eq. (|i~0"D provides the result for the semileptonic decay width b — > cli/i, when 
expressed through the pole b and c quark masses, valid to 0(a 2 ). One notices that Eq. ( |T0"1) 
contains a large second order correction due to the BLM effects. For a long time this obser- 
vation seemed to seriously limit the precision in \ V c b\, as obtained from the inclusive method. 
It was, however, noticed [[1(J that the large value of the second order BLM corrections is 
correlated with the fact that the pole quark masses were used in Eq. (|10D . It is well estab- 
lished, that the pole quark masses can not be defined when non-perturbative corrections 
are addressed. A hint that this is really the case is given by the bad behavior of the pertur- 
bation series itself. In the case of the semileptonic 6-decays the problem is enhanced, since 
the decay width is proportional to the fifth power of the 6-quark mass. 

It was suggested in |16| and then further elaborated in |l5j that the most appropriate 



masses, to be used in the expression for the decay width, are the so-called low-scale running 



1 Actually, ref. || gives — 15.1(a s /7r) 2 . We modify that value by using 4, rather than 3, massless 
flavors for the a s evolution, and by using ^/mbm c , rather than m&, as its normalization point. 
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quark masses normalized at a scale // ~ 1 — 2 GeV. On one hand, such masses can be defined 
on the non-perturbative level, on the other (and this is related to the first point) their use is 
expected to minimize the perturbative corrections to the semileptonic decay width. To 0{a 2 s ) 
a (perturbative) relation between the pole and the low-scale quark masses was obtained in 
11711 and reads: 



m pole = m{ji) + [A(»]p ert + 
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It is important that, in contrast to the pole masses, the low-scale masses do not have 
any significant numerical ambiguity. We use Eq. ([TT|) to rewrite the expression for the 
semileptonic decay width fliTf ) through the low-scale masses normalized at fi — 1 GeV. In 
the BLM approximation such calculation was performed in [[U|. As a result, we find that 
the perturbative coefficients decrease noticeably: 
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where we have used the values of the low scale running quark masses mj = 4.64(5) GeV and 
m c = 1.25(10) GeV, as suggested in [[[J. Though in |TJ the errors assigned to the b and c 
quark low scale running mass were considered conservative, in our opinion this issue is not 
completely clear and a dedicated analysis is required. It is, however, rather certain that in 
contrast to the pole mass, the accuracy of mg can in principle be reliably estimated. 

We see that the perturbative series for the inclusive width appears to behave better when 
the low scale masses are used, in accordance with the theoretical arguments ]16],|T5|,|TJ . 

To sum, we have estimated the second order QCD corrections to the width of the semilep- 
tonic b — > c decay. The small value of these corrections shows that the perturbative series 
is not likely to cause any significant uncertainty in the \V c b\ extracted using the inclusive 
method, provided that the decay width is calculated using the low scale mass definition. 
Further improvement in the theoretical predictions for T s i(B — > X c lv{) will be possible when 
more precise quark mass values and the non-perturbative matrix elements have been deter- 
mined. 

We are indebted to N.G. Uraltsev for many helpful comments. We thank P. Baikov for 
valuable discussions concerning the solution of the recurrence relations for this problem, and 
J.B. Tausk for helping us check some Feynman integrals. This work was supported in part 
by the DOE contract DE-AC02-98CH10886 and by the grant BMBF-057KA92P. 
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FIGURES 




FIG. 1. Status of the two-loop QCD corrections to the decay b — * c + leptons. The dashed line 
denotes the physical region for the actual c quark. Points where the full corrections are known are 
circled. An analytical formula is known along the whole zero recoil line |5],|j . The other two points 
are found from expansions from the base points of the two arrows: at maximal recoil [11] and at 
the intersection with the diagonal (present work). 



